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Abstract 



By exploiting the particle-hole symmetries of the Hubbard model, the 
periodic Anderson model and the Kondo lattice model at half-filling and ap- 
plying a generalized version of Lieb's spin-reflection positivity method, we 
show that the charged gaps of these models are always larger than their spin 
excitation gaps. This theorem confirms the previous results derived by either 

the variational approach or the density renormalization group approach. 
71.27.+a, 71.10.Fd, 75.10.Lp 
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Since the discovery of high temperature superconductivity in the rare-earthy-based cop- 
per oxidesll, interest in the itinerant strongly correlated electron systems has exploded. The 
main concern of physicists is the interplay between the itinerant magnetism and the metallic 
behavior in these systems. As the typical models of the strongly correlated electron systems, 
the Hubbard model!, the periodic Anderson model! and the Kondo lattice mo defl have been 
widely studied in the past several decades. 

To understand the quantum transport and the magnetic properties of the strongly corre- 
lated electron systems, many researcher's interests have focused on the charged gaps and the 
spin excitation gaps in these models at some specific fillings, in particular, at half-filling. For 
instance, by solving exactly the one- dimensional Hubbard model, Lieb and Wu showed that 
the model has a nonvanishing charged gap at half-filling for any on-site Coulomb repulsion 
U > oi. Consequently, in one dimension, the half-filled Hubbard model is an insulator. This 
conclusion is also believed to be true in two dimensions due to the existence of the spin-wave 
excitations, which is caused by the nesting Fermi surface at half-filling, in the system. On 
the other hand, the spin excitation gap of the model is closedil. 

For the periodic Anderson model, the situation is more complicated!,!,!,^. For instance, 
by applying the mean field slave-boson theory, Moller and Wolfle showed! that there exists 
a critical value J c = 4V 2 /U ~ 0.036H 7 , where W is the bandwidth and equal to 4dt for an d- 
dimensional simple cubic lattice, such that, when the hybridization energy V and the on-site 
Coulomb repulsion between <i-electrons U satisfy 4V 2 /U < J c , the periodic Anderson model 
at half-filling has an antiferromagnetic long-range order and the spin excitation gap vanishes. 
On the other hand, when 4V 2 /U > J c , the model has a nonvanishing spin excitation gap 
and its ground state is paramagnetic. The similar conclusions have also been reached by 
numerical calculations^ However, for the Kondo lattice model, such a transition point is 
By using the density-matrix renormalization group method, Yu and White 
found that, in one dimension, the spin excitation gap is nonvanishing for any finite exchange 
interaction J > when the model is half-filledS Their result was further confirmed by 
Shibata et. a/H. 



An interesting observation made by these authors is that: At half-filling, the charged 
gaps of these models are always larger than their spin excitation gapM$^. Therefore, an 
interesting question arose is whether this observation can be re-established on a more rig- 
orous basis. In this paper, by using a generalized version of Lieb's spin-reflection positivity 
technique^PPP, we would like to prove this fact in a mathematically rigorous way As a 
by-product of our proof, one can easily see that this observation is a result of the particle-hole 
symmetry enjoyed by these models at half-filling. 

To begin with, we would first like to introduce several definitions and some useful nota- 
tion. 

Take a finite d- dimensional simple cubic lattice A with an even number of lattice sites 
and impose the periodic boundary condition on it. Then, the Hamiltonian of the Hubbard 
model can be written as 



ieA 



CT <ij> 

+ Uj2(n^-^(n il -^- f iN (1) 

where c[ a {c- i(7 ) is the fermion creation (annihilation) operator which creates (annihilates) 
a fermion with spin a at lattice site i. n- lcr = c\ a Ci a . < ij > denotes a pair of nearest- 
neighbor sites of A. t > and U > are two parameters representing the kinetic energy and 
the screened on-site Coulomb repulsion between fermions, respectively. \l is the chemical 
potential. With respect to Hamiltonian (|l|), the simple cubic lattice is bipartite. In other 
words, it can be split into two separate sublattices A and B such that, the fermion hopping 
takes place only between sites of different sublattices. In the following, we shall exploit this 
fact. 

Similarly, the Hamiltonians of the symmetric periodic Anderson model and the Kondo 
lattice model are given by 

ff A = -iEE ( c tc j(J + c\ a c 

a <ij> 
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In Eq. (|2|), c- la and <ii CT represent the atomic s-orbital and (i-orbital fermion operators at lattice 
site i, respectively. V stands for the hybridization energy of s-electrons and ci-electrons. In 
Eq. (^), o"i and Si represent the spin operators of the itinerant electrons and the localized 
electrons, respectively, and J > is the antiferromagnetic exchange interaction between 
them. 

In terms of Hamiltonians (|2|) and (|3]), the simple cubic lattice is also bipartite. This 
fact can be easily understood by introducing a "double layer lattice structure"El,lii. For 
definiteness, let us consider the periodic Anderson model defined on a specific lattice: The 
two dimensional square lattice with the lattice constant being set to be unit. We take two 
identical copies of this lattice, Ax and A 2 , and make a doubly-layered lattice A by connecting 
the corresponding lattice points of Ai and A 2 with bonds of length a = 1. Now, each point 
of A has coordinates r = (i, m) with m — 1, 2. Obviously, A has 2Aa lattice points. Next, 
we define new fermion operators e ra by 

c ia , if m = 1; 

(4) 

d- 1(T1 if m = 2. 

With the definitions of A and e rCT , the Hamiltonian of the periodic Anderson model can 
be thought as the Hamiltonian of a generalized Hubbard model on the bipartite lattice A, 
if V is taken to be the "hopping energy" of e-electrons between layer 1 and layer 2. 

Similarly, for the Kondo lattice, if we re-define the partition of sublattices by the following 
rule: The hopping energy t and the exchange energy J may be nonzero only for a pair of 
sites belonging, respectively, to the different sublattices, then lattice A as well as the original 
lattice A are apparently bipartite in terms of Hamiltonian (pi) . 



The Hamiltonians Hh, Ha and Hk enjoy several symmetries, which make the analysis 
of these models easier. 

First, Hamiltonians ([!]), (0) and (|3|) commute with their total particle number operators 
N, respectively Consequently, their Hilbert spaces can be divided into numerous subspaces 
{V(N)}. Each of these subspaces is characterized by an integer N, the total number of 
fermions in the system. In particular, the subspace V(N) is called the half-filled subspace 
if N — N\ for the Hubbard model, iV = 2N\ for both the periodic Anderson model and the 
Kondo lattice model. 

Furthermore, if we define the total spin operators, for the Hubbard model, by 

= \ E ( C iT C U + C U C it) ' 
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for the periodic Anderson model, by 

4 = ^E ( c it c u + c u c iT + d h d n + 4i<h) , 

z ieA 

4/ = 97 Z ( c iT c U - c i| c iT + d h d ii - dl^) , 
ZL ieA 

S, = ~£(>n (6) 
z ieA 

and, for the Kondo lattice model, by 

Sx = \y. (4,^ + cj |CiT + fy n + 4/ iT ) , 

z ieA 

Sy = 97 Z ( c lt c u - c u c n + /it/u - /i/n) > 

Zi ieA 

& = 9E("iT-"U + ^T- n u) ( 7 ) 
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then it is easy to check that the Hamiltonians H^, Ha and Hx commute with their cor- 
responding total spin operators S + = S x + iS y , S- = S x — iS y and S z . Consequently, any 
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eigenstate \l/ n with quantum number S 2 = S(S + 1) of these Hamiltonians must have 2S + 1 
isotopes {* n (M)} with —S < M < S. 

Here, we would like to emphasize that, in Eq. (|7|), fl a and f- 1(T represent the fermion 
operators of the localized spins. Consequently, they should satisfy the following constraint 
conditionil 

AVAt + flih = i (8) 

This makes them different from operators c ia and c- ia , which are the itinerant electron oper- 
ators. 

When the systems are half-filled, the chemical potential /i = for all the Hamiltonians 
Hh, Ha and Hk- This fact is due to the particle-hole symmetry enjoyed by these models 
at the special fillingiil. As a result, the Hamiltonians Hh, Ha and Hk also commute with 
the so-called pseudo-spin operators, which are defined, for the Hubbard mo del0, by 

J+ = E e (i) c iT c u, J- = J+, 
ieA 

^=«E("it+%- 1 ) ( 9 ) 

z ieA 

for the symmetric periodic Anderson model^, by 

ieA 

^^EK + ™ii + < + <- 2 ) ( 1Q ) 

z ieA 

and, for the Kondo lattice mo del@, by 

i+ = Ee(i)(W-/iVi). J- = Jl, 

ieA 

J* = o E (<t + n u + 4 + n d ~ 2 ) ( n ) 

z ieA 

where e(i) = 1, fori £ A; e(i) = —1, fori £ B, is the alternating function. These operators 
also satisfy the commutation relations of the conventional spin operators. Since Hh, Ha 
and Hk commute with J + , J_ and J z , both J 2 and J z are also good quantum numbers and 
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hence, each eigenstate of Hh-, Ha and Hk is characterized by a quantum number J and a 
quantum number J z with —J<J Z <J. It has been shown that, at half- filling, the ground 
states of the Hubbard mode$, the symmetric periodic Anderson rnode^ and the Kondo 
lattice mo<ie@ ; @ on the simple cubic lattice have quantum numbers 5 = and J = 0. 

Next, we would like to introduce the definitions of the charged gaps A c and the spin 
excitation gaps A s for these strongly-correlated electron models at half-filling. In Refs.@ 
and§, these quantities are defined by 

A c = E (J = 1, S = 0) - E (J = 0, S = 0), 

A s = E (J = 0, S = 1) - E (J = 0, S = 0) (12) 

where Eo(J = j, S = s) is the lowest eigenvalue of the corresponding Hamiltonian in the 
subspace with quantum numbers J = j and S = s. By using the definitions of the pseudo- 
spin operators and considering the fact that the ground states of these Hamiltonians at 
half-filling have quantum numbers S = and J = 0, the above definitions can be also 
re-written into the following formsEl, in terms of the total fermion numbers 

A c ee E (N + 1) + E (N - 1) - 2E (N), 

A s ee Eq{N, S = l)- Eq(N, S = 0) (13) 

where N = N\ for the Hubbard model and iV = 2N\ for both the periodic Anderson model 
and the Kondo lattice model. Notice that, at half-filling, the Hamiltonians Hh, Ha and 
H K enjoy the particle- hole symmetry and hence, identity E (N + 1) = E (N -1) holdsS. 
Therefore, A r can be further written as 



A c = 2 



E (N + 1)-E (N) (14) 



With these preparations, we shall now summarize our main result in the following theo- 
rem. 

Theorem: For the Hamiltonians Hh, Ha and Hk defined on a rf-dimensional simple 
cubic lattice, when the system is half-filled, the charged gaps and the corresponding spin 
excitation gaps satisfy the following inequality 



A c > A s (15) 

Proof of the theorem: To prove this theorem, we shall apply a generalized version of 
Lieb's spin-reflection positivity methodli,0, which we previously used to study the binding 
energy of fermions in the negative- U Hubbard mo deM In order to make it more readable, 
we shall divide the proof of this theorem in several steps: 

(1) First, by introducing a unitary partial particle-hole transformation for each model, 
we map the original Hamiltonians to some equivalent Hamiltonians with negative coupling 
constants. 

(2) Then, we write each of the transformed Hamiltonians into a form of the direct product 
of up-spin fermion operators and down-spin fermion operators. 

(3) We shall then apply the spin-reflection positivity method to these transformed Hamil- 
tonians and establish an inequality for the lowest eigenvalues of these transformed Hamilto- 
nians in the different subspaces. 

(4) Finally, we apply the inverse of the partial particle-hole transformations to the in- 
equality and finish the proof of the theorem. 

In the following, to avoid unnecessary digression, we may directly apply some well- 
established mathematical results without proving them. Naturally, in that case, we shall 
refer to some standard references for the reader's convenience. 

We now proceed to the proof of the theorem. 

Step 1: It is a well-known fact that, for the Hamiltonians H^, Ha and Hx at half-filling, 
there exist unitary transformations Uh, Ua and Uk, which are called the partial particle- hole 
transformations^. Under these transformations, each Hamiltonian with positive interaction 
coupling constants is mapped to a corresponding Hamiltonian with negative interactions. 
To be more precise, let us consider these Hamiltonians one by one. 

(a) For the Hubbard Hamiltonian ([]]), U H is defined 

U 1 h cvUh = At, U ] „c u U„ = £ (i)4 (16) 
At half-filling (/i = 0), under U H , H H is mapped to 



H H = uIh h u b = -tJ2 E ( c U. + 4^< 

CT <ij> 

^g(»u-i)(«u-i) (17) 

Apparently, if// has the same form as Hfj. However, in Hamiltonian (|17|) , the sign of £/ is 
changed. 

(b) For the symmetric periodic Anderson model, the unitary transformation U a is defined 
by@ 

U\c^U A = c iT , UaChUa = e(i)4|, 

U\f- A U A = /i T , UaM&a = -e(i)/J (18) 
Consequently, under U A , Ha is mapped to 

H A = U a H a Ua = -t E E + 4 a c ia ) 

CT <ij> 

+ ^EE( c U. + 4c i( x) 

Notice that the sign of U is changed but the sign of V is invariant. Since the hybridization 
term can be mathematically treated as a generalized "hopping" term, as we shall show, the 
sign of V plays no role in the following proof. 

(c) For the Kondo lattice model, the unitary partial particle-hole transformation Uk is 
given by0 

U K c- A U K = d T , U K c n U K = e(i)4j., 

Uy^U K = f ih U K f n U K = -e(i)fl (20) 
Under Uk, the constraint condition (||) now reads 

./r.A- = (21) 
and the half-filled Hamiltonian Hk is mapped tc 



H K = UlH K U K = -tJ2 E ( C U- + 4*°*") 

u <ij> 

— T E E ( C icr C io- + fiafiv) 
4 cr ieA 

+ 9 E ( C iT C it/n/iT + C i| C ii/ii/i|) 

z ieA 

-^EKM/u+zMiPu) (22) 

z ieA 

We notice that, in Hamiltonian (P^)> the sign of the last term is negative. This fact is the 
basis of our proof of the theorem for the Kondo lattice model. The first three terms can 
be mathematically treated as generalized "hopping" terms. We shall see that their signs do 
not matter. 

Here, we would like to make some remarks. 

Remark 1: Under U H , Ua and Uk, the half-filled subspace for the corresponding Hamil- 
tonian is mapped into itself. In particular, since these transformations are unitary, the 
ground states of Hh, Ha and Hk in the half-filled subspace V(N) are mapped onto their 
counterparts for H H , Ha and H K in the same subspace. However, other subspaces are not 
invariant under these transformations. 

Remark 2: Under the partial particle-hole transformations, the spin operators S + , 
S- and S z related to each Hamiltonian are mapped onto the corresponding pseudo-spin 
operators J + , J_ and J z , and vice versa. Consequently, under these transformations, an 
eigenstate J = j, S = s) of the original Hamiltonians H H A , k is mapped onto an eigenstate 
^(J — s, S — j) of the transformed Hamiltonians Hh,a,k- In particular, the ground states 
of the original Hamiltonians in the sector with quantum numbers J = j, S = s is mapped 
to the ground states of the transformed Hamiltonians in the sector with quantum numbers 
J = s, S = j. This observation is important in carrying out our proof in Step 4. 

Step 2. Next, we would like to write Hamiltonians (0); (0) an d ( p2|) into a form of 
the direct product of up-spin fermion operators with down-spin fermion operators. For this 
purpose, we shall follow Ref.0 and introduce the following new fermion operators for each 
Hamiltonian. We let 
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A T = e iT , C n = (-l) JV t eu (23) 



where ei CT stands for c^, and c? io . appearing in Eqs. (|17), ( |19|) and (p2|) , respectively. 
iVf is the total number of fermions with up-spin in the system. Here, we would like to 
emphasize that the new fermion operators {Ci|}j now, commute with {At}- Consequently, 
Hamiltonians (|T7|) , (|T9| ) and (|22|) can be, respectively, re- written as 

-^E^-")® (n U -~) (24) 



H A = f'®i + / ® f[ 

-^E (25) 



and 



E K = T'{ ® J + J ® Tf 

- | E (4 T /iT ® cl^ + // T c iT ® /.\c u ) (26) 
z ieA 

In Eqs. (0), (p5|) and fl2"6|) , T CT , and T" are some Hermitian polynomials of the fermion 
operators of spin a. I is the identity operator. Apparently, all the above Hamiltonians can 
be written in the following standard form 

H = G^®i + I®G i 

-A£(&T®Qii + $ t ®Qu) (27) 
ieA 

where A is a positive constant. We notice that all the operators {Qi a } and {<0L-} i n Eq. fl2"7|) 
are real operators. In other words, they are polynomials of C\ a and C\ a with real coefficients. 
The fact is of fundamental importance for applying Lieb's spin-reflection positivity method. 
Now, we are able to treat Hh, Ha and Hk simultaneously, by studying the standard Hamil- 
tonian 
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Step 3: Now, let us consider the ground state \I/o(jV + 1) of H in the subspace V(N + 1). 
Since the spin operators S + and S- commute with H, by applying these operators an 
appropriate number of times, we can always transform \l/ (^+ 1) into a state satisfying the 
condition jV-j- — iVj_ = 1. This state has quantum number S z = ~. In the following, we shall 
exclusively use *$?o(N + 1) to denote this state. 

The wave function ^o(N + 1), which has y + 1 up-spin fermions and ~ down-spin 
fermions, can be written as 

y (N+l) = Y,W mnX L®X l n (28) 

m, n 

In Eq. (|28|), x% is a state vector defined by 

xl = C} ia ---C} Ma \0) (29) 

where (ij, . . . , ij^f ), M = y + 1, for a =|; M = y, for cr =|, denote the positions of fermions 
with spin cr. We would like to emphasize that, when H = H^, the constraint condition 



should be taken into consideration to determine the coefficients {W^mn} of the ground 
state wave function. Apparently, the entire set {xl} gives a natural basis for V a (M), the 
subspace of M fermions with spin cr. However, we should notice that, if we naively choose 
TCi = V|(y + 1) and Tii = Vj_(y) to be the subspaces for up-spin and down-spin fermions, 
then the coefficient matrix W = (W mn ) will be an C^ 2+l x C^ 2 matrix, which is not a 
square matrix. Mathematically, it is rather difficult to deal with such a matrix. To avoid 
this nuisance, we shall define both Ti-\ and Hi by H a = K-(y) ©V^(y + 1). Consequently, 
Til an d have the same dimension and hence, matrix W can now be written as an D x D 
square matrix with D = C^ 2 + C~^ 2+1 . Explicitly, we have 



W 



( O 

\° I 



(30) 



In Eq. (|30|) , M. is an C^ 2+1 x C~^ 2 nonzero matrix and O represents zero matrices. In 
terms of matrix W, the normalization of ^q(N + 1) is now given by TrW^W = 1. 
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For such a square matrix, we have the following polar factorization lemma in matrix 
theory. 

Lemma: Let A be an arbitrary (Not necessarily Hermitian) n x n matrix. Then, there 
are two n x n unitary matrices U, V and annxn diagonal semi-positive definite matrix H 
such that 

A = UHV, 

h mn = h m 5 mn and h m > 0, m = 1, • • • , n. (31) 

One can find the proof of this lemma in a standard textbook of matrix theory^ or read 
the appendix of Ref.0. 

Applying this lemma, we can find two unitary matrices U, V and a diagonal positive 
semidefinite matrix H, such that W = UHV. Consequently, ^o(N + 1) can be re-written 



as 



%(N + 1) = E E W mnXl ® xl 
m=l n=l 

= EE (uHV)mn xl ® xi 

m=l n=l 

= f>V>J®0} (32) 
i=i 



with 

D D 



= E ^xL ^ = E ^ (33) 



m=l n=l 



Since [/ and V are unitary, } and {4>j} are also orthonormal bases in subspaces Ti^ and 
Tii, respectively. Furthermore, since ^o(N + 1) is an eigenvector of iVj, we have 



iV T | ^ (iV + 1)) = + l) | ^ (iV + 1)) (34) 



or, equivalently, 



D 

I 



E hi [aW 

1=1 

D 

Eh 

i=i 
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i (35) 



Taking the inner product of Eq. fl35|) with <p\, projects out 



hp 



N 
~2 



(36) 



Consequently, the corresponding state if)J is an eigenvector of iVj with eigenvalue N/2 + 1, 
if hi 7^ 0. The same conclusion can also be reached for operator iVj_ and states {0/"}. 

For technical reasons, the above conclusions are generally written in the following weaker 
form 



1) | iV T | V (N+1)) 
D ~ N 

i=i A 



(37) 



and 



(tf„(W+l) I I ^o(A^+l)} 



iV 



(38) 



In both Eqs. fl37|) and (0), the spin indices are dropped in the sums, because, in each 
equation, only one species of spin is involved. These equations are particularly useful in 
proving the strictness of inequality ([T^) for a finite system (See Ref.0 for details) and will 
be used in the following. 

In terms of Eq. (|3~2|) , the ground state energy of H in subspace V(N + 1) is now given 

by 

E (N + 1) 
= (* (N + 1) I H | V (N + 1)) 
= £> a \(if>J | G T I if>}) + (4>\ I Gi | <j>\) 



i=i 



D 



aE E I Git 1 4w I 

ieA Vx, i 2 =l 
ieA Vi,z 2 =i 



4) 



0i) 



(39) 
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Applying inequality |a6|<|(|a| 2 +|6| 2 )to each term in the triple summations and 
dropping the spin indices, we obtain 

E (N+1) 

>^X>* \(ipi\G\ipi) + (ipi\G\ipi) 
A 1=1 

+ \j2 h i [(h \ d \ + (<h \ & \ h) 
1 1=1 

D 



^E( E h 1 h h (ijj h \Qi\ilJ h )(il) h \Qi\'il)i 1 ) 
z ieA \h,i 2 =i 



D 



IE! E Khh{<t>h I Qi I <t>h){<l>h I Qi I <t>h) 

z ieA \«i,«2=i 



D 



^E( E KhMiAQ\\^)^iAQ\\^) 

z ieA \«i,«2=i 



D 



/ ieA \«i,«2=i 

Next, we introduce new wave functions ^\ and \J/ 2 by 

D D 

*i = E M>/ ® ^ 2 = E 

!=1 i=l 



(40) 



(41) 



where and 0; are the complex conjugates of and 0/, respectively. Apparently, we have 



]T/i 2 = (^o(iV+l)|^o(iV + l)> = l 



(42) 



Since G is hermitian and {Qi} ({Q\}) are real, in terms of and ^2, inequality 
be re-written as 



can 



E (N + 1) > \H\^ 1 ) + ^ 2 \H\ * 2 > 



(43) 



On the other hand, we notice that ^1 and ^2 are actually wave functions in subspace 
V(y + 1, y + 1) and V(y, y), respectively. For example, by using the constraint condition 
(1371), we have 
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D iV 



1=1 / 
Therefore, by the variational principle, we obtain 

E (N+ 1) > ^E Q (N + 2) + l -E a (N) (45) 

Remark 3: We should notice that, while the construction of and ^ for both Hh 
and Ha is straightforward, the job should be done in a more cautious way as far as Hk is 
concerned. We should show that both and ^2 satisfy the constraint condition (|21|) . That 
can be achieved in the following way 

By reorganizing the rows and columns of the coefficient matrix W (by a unitary transfor- 
mation), we can always write it in a "block diagonal" form: W = diag (Wi, W2, • • •)> where 
each block is a square matrix and corresponds to a sector in which the distribution of 
/-fermions is specified, subject to the condition fl2~l~|) . Naturally, Wk has also the form of 
Eq. (|30l). Applying the lemma to each of these submatrices {V\4}, we can construct unitary 
matrices U and V with the following block diagonal form: 

U = diag (U u U 2 ,- ■ •) , V = diag (Vi, V 2 , ■ ■ •) (46) 

and diagonalize the coefficient matrix W by them. Now, it is easy to see that the wave 
functions and ^2 constructed according to the above-mentioned rule satisfy the constraint 



condition (|2lD . 

Step 4-' Finally, we apply the inverse of Uh, Ua arid Uk to map the "negative coupling" 
Hamiltonians H^, Ha and Hx back onto the original Hamiltonians Hh, Ha and Hk- We 
are mainly interested in the change of inequality (pE5[) under these transformations. For this 
purpose, let us consider how each of E (N), E (N + 1) and E (N + 2) changes under the 
unitary inverse transformations. 

First, as we mentioned in Remark 2, Eq(N) is unchanged under the partial particle- 
hole transformations. In other words, it still presents the lowest eigenvalues of the original 
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Hamiltonians Hh, Ha and Hk in the half-filled subspace V(N = N). Interestingly, Eq(N+1) 
is also invariant. This is due to the fact that the partial particle-hole transformations only 
change the particle number of down-spin fermions from | to Af - | = | and keep the 
particle number of up-spin fermions, = y + 1 unchanged. Consequently, the subspace 
V(iV| = — + 1, iVj. = —) is unitarily mapped into itself. However, the change of E (N + 2) 
demands a careful consideration. 

It has been previously proven that, in the subspace V(N + 2), the ground states of 
the Hamiltonians iijlljill, ifjil and ifft-@,IH have quantum numbers S — and J — I. 
Therefore, under the inverse transformations U^ 1 , U^ 1 and U^ 1 , these states are mapped 
onto the ground states of the original Hamiltonians Hh, Ha and Hk in the subspace with 
quantum numbers J = and S = 1, as we discussed above in Remark 2. In other words, 
these states are the ground states of the original Hamiltonians in the half-filled sector with 
S — 1. Consequently, we obtain 

E (N + 2; H h ,a,k) = E (N } 3 = 1; H H>AjK ) (47) 

Therefore, for the Hamiltonians H H , Ha and H K , inequality fl4"5p now reads 

E (N + 1) > ^E (N, 5 = 1) + ~£ (AT, 5 = 0) (48) 

or, equivalently 

E (N + 1)-E (N) 
> ^E (N, 5 = 1)- ±Eo(N, 5 = 0) (49) 

Multiplying both sides of inequality fl49|) by 2, we obtain inequality (p!5[). That ends our 
proof of the theorem. QED. 
Some remarks are in order. 

Remark 4: Actually, inequality ( [49"D can be slightly strengthened for the Hubbard model 
and the symmetric periodic Anderson model at half-filling. Following Ref.0, one can easily 
show that, for any finite lattice A, inequality (f49D is strict. In other words, the systems have 
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a Mott charged gap. However, in the thermodynamic limit, the inequality resumes the form 
of Eq. P§p . Therefore, to study the Mott metal-insulator transition of these models in the 
thermodynamic limit, one needs to introduce some more sophisticated method. 

Remark 5: In the above proof, we assumed that the coupling constants in each Hamil- 
tonian are site- independent. That is only a harmless assumption to make our proof simpler. 
Actually, inequalities (|45|) as well as ([49]) still hold for the models with site- dependent cou- 
pling constants and its proof is essentially unchanged. One can find more details in Ref.0. 

In summary, in this article, by exploiting the partial particle- hole symmetry of the Hub- 
bard model, the periodic Anderson model and the Kondo lattice model at half-filling and 
applying a generalized version of Lieb's spin-reflection positivity method, we proved that 
the charged gaps of these models are always larger than their spin excitation gaps. This the- 
orem confirms the previous results derived by either the variational approach or the density 
renormalization group approach. 
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